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Abstract: The focus point region of super symmetric models is compelling in that it simul- 
taneously features low fine-tuning, provides a decoupling solution to the SUSY flavor and 
CP problems, suppresses proton decay rates and can accommodate the WMAP measured 
cold dark matter (DM) relic density through a mixed bino-higgsino dark matter particle. 
We present the focus point region in terms of a weak scale parameterization, which al- 
lows for a relatively model independent compilation of phenomenological constraints and 
prospects. We present direct and indirect neutralino dark matter detection rates for two 
different halo density profiles, and show that prospects for direct DM detection and indirect 
detection via neutrino telescopes such as IceCube and anti-deuteron searches by GAPS are 
especially promising. We also present LHC reach prospects via gluino and squark cascade 
decay searches, and also via clean trilepton signatures arising from chargino-neutralino pro- 
duction. Both methods provide a reach out to rUg ~ 1.7 TeV. At a TeV-scale linear e~^e~ 
collider (LC), the maximal reach is attained in the Z1Z2 or channels. In the DM 

allowed region of parameter space, a ^/s = 0.5 TeV LC has a reach which is comparable to 
that of the LHC. However, the reach of a 1 TeV LC extends out to rrig ^ 3.5 TeV. 
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1. Introduction 

Supersymmetric models of particle physics provide a compelling case for physics beyond the 
Standard Model (SM). However, in spite of their many successes, they also provide a long 
list of potential problems. For instance, supersymmetric models are supposed to provide a 
solution to the fine-tuning problem which arises when the SM is embedded in theories with 
high mass scales beyond a few TeV. However, the generally excellent agreement of precision 
EW observables with SM predictions, along with null search results for various rare decays 
and other loop induced processes, points to a rather heavy sparticle mass spectrum, with 
sparticles typically in the TeV regime. These observations are supported by recent search 
results from LEP2 that the chargino mass > 103.5 GeV and the SM Higgs mass 

rriHsM > 114.4 GeV. The latter limit, when applied to the Higgs bosons of the MSSM, also 
implies relatively heavy top squarks. A rather heavy sparticle mass spectrum, on the other 
hand, naively seems to re-introduce the fine-tuning problem into supersymmetric models. 
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In addition, in the 124 parameter Minimal Super symmetric Standard Model (MSSM), 
unsuppressed FCNC effects arise from general Lagrangian parameters, as do large contri- 
butions to the electric dipole moments of the electron and neutron from possibly large CP- 
violating phases. It has been noted by many authors that scalar masses in the multi-TeV 
regime can act to suppress most or all of these undesired effects to within levels tolerated 
by data]^]. Multi-TeV scalar masses can be reconciled with a no-fine-tuning requirement in 
two cases. In one case, inverted scalar mass hierarchy models|Q] (IMH) require scalars of the 
first two generations to be in the multi-TeV regime, while scalars of the third generation, 
which enter fine-tuning calculations, remain at sub-TeV levels. In practice, the radiatively 
driven IMH models require t — b — t Yukawa coupling unification and non-universal soft 
SUSY breaking (SSB) Higgs masses to be viablej^, |^. The second case, the topic of this 
paper, is that of hyperbolic branch||5| or focus point modelsQ (HB/FP), wherein all three 
generations of scalars can be in the multi-TeV regime, while fine-tuning is respected for 
low values of the GUT scale universal gaugino mass 772-1/2- 

The HB/FP region appears already in Ref. Q as a region in the tuq vs. mi/2 plane 
of the minimal supergravity (mSUGRA) model where mo is in the multi-TeV regime, but 
where the absolute value of the superpotential /x parameter becomes small, adjacent to 
regions where radiative electroweak symmetry breaking (REWSB) fails to occur (where 
jj? < 0). The small value of |^| leads to a mixed higgsino-bino LSP (a mixed higgsino dark 
matter (DM) candidate) and a rather light spectrum of charginos and neutralinos. This 
mixed higgsino region with large scalar masses was investigated more thoroughly by Chan, 
Chattopadhyay and Nath in Ref. |^], where it was noted that the ^ parameter can be 
regarded as a measure of fine-tuning, and where the trajectory of constant /x was found to 
form a hyperbolic branch trajectory in the mSUGRA parameter space. In Ref. [^], Feng, 
Matchev and Moroi noted the focusing behavior of the rn'jj^ renormalization group (RG) 
trajectory, wherein a variety of GUT scale iTi'jj^ values would be "focused" to a common 
weak scale value of 'm'jj^. Since REWSB leads to /i^ ~ —mj^^ at the weak scale, the focused 
solutions gave rise to small fj? values. These authors moreover performed a sophisticated 
fine-tuning analysis, and showed that fine-tuning was small in the HB/FP region as long as 
mi/2 was not too large. It could be seen in Ref. and more fully in |^, IC, 0], that the 



relic density is indeed low in the HB/FP region, as is typical of mixed higgsino dark matter. 



Further, in Ref. Q, it was shown that a variety of directjl^, 13| and indirect[14, 15| DM 
detection rates were large in this region, due to the large neutralino-nucleon scattering cross 
sections, and also due to the large neutralino-neutralino annihilation rates. In addition, 
collider reaches in the HB/FP region were found for the Fermilab Tevatron[p!^], the CERN 
LHC0, 0, |18| and the International Linear Colhder (ILC)|l|]. 



All of the above HB/FP analysis were performed in the niQ vs. mi/2 plane of the 



mSUGRA model. However, it was noted in Ref. [16| that the exact location of the HB/FP 
region in the mSUGRA plane is extremely sensitive to the assumed value of mj . In addition, 
different algorithms for predicting sparticle masses in the mSUGRA model were found to 
give very different portrayals of the shape and location of the HB/FP region even for the 
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same assumed value of mi[20|^. Finally, it was shown in Ref. |21, ^ that the HB/FP 
region occurs also in models with non- universal scalar masses, and that its location can 
shift depending on the amount and type of non-universality which is assumed. 

In this paper, we propose a more sensible parametrization of the HB/FP region based 
purely on weak scale parameters fj, and U{1) gaugino mass Mi. Such a parametrization 
should be independent of the GUT-to-weak scale evolution algorithm assumed, and also 
should be only weakly dependent on the value of mt or other SUGRA parameters which 
are assumed. In Sec. we present our re-parametrization of the HB/FP region, and show 
the regions of the weak scale Mi vs. /i plane which are allowed by the recent WMAP 
measurement ClcDAih'^ = 0.113 it 0.009 of the relic density of cold dark matter in the 
universe. In Sec. |3|, we show current astrophysical constraints on the HB/FP region 
arising from various sources, including an overproduction of ^Li in the early universe. 
We also show the prospects for exploring the HB/FP region via direct and indirect DM 
detection. In the HB/FP region, prospects are especially encouraging via Stage 3 direct 
detection experiments, detection of DM via neutrino telescopes, and via anti-deuteron 
searches at experiments such as GAPS. The latter results depend strongly on the galactic 
DM density profile which is assumed. In Sec. ^, we show the reach prospects of the CERN 
LHC and also for the ILC for SUSY in the WMAP allowed part of the HB/FP region. The 
CERN LHC reach is evaluated in the case of gluino and squark cascade decays, and also 
for clean trileptons arising from W1Z2 production. In both cases, the LHC reach is out to 
rrig ~ 1.7 TeV assuming 100 fb^^ of integrated luminosity. The ILC reach is maximal in 
the e~^e~ Z1Z2 or Z1Z3 channels, and extends out to nig ~ 1.7 TeV for a ^/s = 0.5 TeV 
ILC, but to rUg ~ 3.5 TeV for a ^ = 1 TeV colhder. 

2. Unraveling the Focus Point Region 

The purpose of this section is to introduce the HB/FP parameter space region of the 
mSUGRA model which produces a sufficiently low relic neutralino abundance in the Early 
Universe at large values of the universal scalar mass mo- We motivate here why a GUT- 
scale description of the HB/FP region is ill-suited for phenomenological studies, and outline 
an alternative and complementary weak scale parameterization. The latter allows us to de- 
scribe in detail the cosmologically allowed areas in the more physical weak-scale parameter 
space, and will be used throughout the remainder of this report to study the corresponding 
phenomenology at dark matter search experiments and at colliders. 

Radiative electroweak symmetry breaking is the mechanism in which EWSB is trig- 
gered by rnjj^ turning negative due to its RG evolution. The RGE for rnjj^ reads 

where 



Xt = (m^3 + ml^ + ml^ + A^) , (2.2) 

S = ml,, - ml^ + Tr [m^ - mi - 2ml + + m|] , (2.3) 

^Even updated versions of the same computer codes would give shifted locations of the HB/FP region. 
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where 5 = at all scales in models with universality. Among the necessary conditions for 
the spontaneous breaking of the EWS, the one which determines the (absolute) value of 
the n parameter reads (at tree level) 

^ ~ tan2/5-l 2 ' ^^'^^ 

At moderate to large values of tan /3, as implied by the LEP2 lower bound on the lightest 
Higgs mass, /i^ ~ —mjj^. When the universal mSUGRA GUT scale scalar mass mo 
takes values much larger than all other soft-breaking masses, the RG evolution of m'jj^ is 
dominated by the term Xt in Eq. (^.Ij). As mo increases, cancellations occurring in Xt yield 
smaller and smaller absolute values for m^^ . Therefore, Eq. ( [2.4]) leads to the possibility of 
achieving, in principle, arbitrarily low values for /i^, until, eventually, < 0, and REWSB 
can no longer be obtained. 

Within the mSUGRA model, low values of /i therefore occur in the region of very 
large universal scalar mass mo ^ 1 TeV. This opens up a qualitatively new window in the 
model parameter space, as far as the cosmological abundance of thermally produced relic 
neutralinos is concerned. As /i approaches the low-scale value of the lightest soft-breaking 
gaugino mass, Mi in the case of mSUGRA, the higgsino component of the LSP increases, 
leading to an enhancement of efficient LSP pair annihilations into gauge bosons (the latter 
being largely suppressed in the case of a bino-like LSP). Further, the lightest chargino and 
the two next-to-lightest neutralinos get closer in mass to the LSP, contributing as well to 
the suppression of the LSP relic abundance through co- annihilations. 

As a result, the particle spectrum of the cosmologically allowed HB/FP region of the 
mSUGRA model is characterized by (1) a heavy scalar sector, in the multi-TeV range (with 
the exception of the lightest CP-even Higgs boson), and (2) low values (anywhere below 
or around 1 TeV) of the n parameter. For ^ ^ 1 TeV, the requirement of a sufficiently 
low neutralino relic density forces moreover the relation ~ Mi, the latter being the 
soft-breaking hypercharge gaugino mass. The weak-scale values of the three soft-breaking 
gaugino masses, unified to a universal value mi/2 at the grand unification (GUT) scale, are 
given by the usual GUT relations. 

Ml ~ 0.44 mi/2 (2.5) 

Af3 = M, (2.7) 

Evidently, the decoupling of the supersymmetric sfermion and heavy Higgs sectors 
implies that the critical parameters entering the phenomenology of the HB/FP region at 
colliders and at dark matter search experiments are Mi and ^, their relative size setting 
the LSP mass and higgsino fraction. Further, the gluino mass is determined as well by Mi, 
through Eqs. (|2.5| )-( ^?7|) . The sign of /i and tan/5 affect as well, though less critically, the 
low-energy implications of the setup. 
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Figure 1: The projection, onto the (Afi,/i) plane, of some mSUGRA parameter space sUces at 
constant mo, for rrit = 174.3 GeV, = 0, /i > 0, at tan/3 = 10 (left panel) and 50 (right panel). 
The blue lines indicate the points featuring the maximal value of at fixed mi/2. 



In the standard GUT-scale parameterization, where one slices the parameter space e.g. 
along (m-o 1^1/2) planes, the HB/FP region appears as a narrow line squeezed on the large 
niQ region adjacent to where no REWSB is attainable. The steep and complicated behavior 
of /X as a function of mo makes it difficult to read out the neutralino mass and composition 
in that representation. Further, the HB/FP region is often plagued by numerical stability 
problems, which will be discussed in detail in Sec. 2.1, affecting the phenomenologically 
crucial low-scale value of the /x parameter. As a consequence, it is not easy to read out 
from the standard (w-q, parameterization most of the phenomenologically relevant 

information, e.g. the neutralino mass range and composition compatible with the WMAP 
relic abundance, and the projected reach of collider and dark matter search experiments. 

The purpose of our analysis is therefore to map the GUT-scale representation of 
mSUGRA in terms of the universal scalar and gaugino masses onto the more physical 
(Mi,/i) plane on which the whole HB/FP region phenomenology sensitively depends. 

From the high-energy scale point of view. Mi depends essentially only on mi/2 (see 
Eq. (p.5|)), while /x is sensitive, in principle, to all mSUGRA GUT input parameters. For 
given fixed values of tan (3 and of the trilinear coupling ^0 > the mSUGRA parameter space 
spans a limited region only, on the physical (Mi,^) plane. The // parameter features, in 
fact, a certain maximal value (reached at low or intermediate values of mo) along slices at 
fixed mi/2- Values larger than the mentioned maximum cannot be obtained if REWSB is 
required. To illustrate this point, we show, in Fig. 0, curves at fixed mo on the (Mi,/i) 
plane, computed with Isajet 7.72 ||2^ at two values of tan /5 = 10 and 50, with an input top 
mass mt = 174.3 GeV. The iso-mo curves terminate where the RG code no longer converges 
to a stable solution, while the blue curves indicate the maximal /i values. We remark that 
the low /i termination point of the curves is in fact only a numerical artifact, and that 
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the entire low // region extending to = is a physically viable portion of the parameter 
space, which is not accessible in many codes which rely on a GUT- scale parameterization. 

As expected from Eq. (p.4|), a given value of ^ is obtained, at smaller tan/3, with 
larger values of mo, and vice- versa: this reflects the well-known fact that the HB/FP 
region appears at smaller values of mo the larger tan (3 is. 

The sensitivity of the HB/FP region parameter space in the (mo, mi/2) plane on 



the value of mt has been widely reported, see e.g. Ref. fig, g4|, |2^, and can be readily 
understood from Eq. (|2.1[) . The shape of the weak scale parameter space (the (Mi,/i) 
plane) will change as well, but is much less sensitive to the assumed value of mt. We show 
the analogue of Fig. |l| but at mt = 178 GeV in Fig. 

2.1 Numerical issues in the HB/FP region 

As we outlined above, a first, and critical, numerical issue in the HB/FP region regards the 
possibility of achieving a stable and convergent solution for low values of the parameter. 
To investigate this problem, and to verify the consistency of numerical results in the pa- 
rameter space region of interest, we compared the latest release of Isajet, v. 7. 72, with the 
latest release of another RGE evolution package. Suspect, v. 2.34 (for a comparison among 
the results of these and other numerical codes, see also Ref. p^, |2j, |20|). 

Fig. ^ illustrates how the two RG evolution codes project iso-mo lines onto the physical 
parameter space of the focus point region, the (Mi,/i) plane. We adopt for both numerical 
codes the same input top mass, which we set to mj = 174.3 GeV. The lines end, as in 
Fig. |l| and ^, where no stable solutions are found, in the low /i portion of the plots. We see 
that in both cases, values of /i smaller than 100- 200 GeV in the large mi/2 region cannot 
be numerically resolved. Further, although both RG codes feature 2-loop RG running of 
gauge and Yukawa couplings, the disagreement is remarkable (particularly in the low ^ 
region and at low values of tan/3). 
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Figure 3: A comparison, along the same parameter space slices as in Fig. y, between the two RG 
evolution nmnerical codes Isajet 7.72 [H (red dashed hnes) and SUSPECT 2.34 fH] (black soHd 
lines) 



The calculated value of the parameter as a function of mo along mSUGRA slices at 
fixed is shown in Fig. ^. While the agreement among the two codes is excellent in the 
low mo end of the plots, when approaching the HB/FP region the values of /U (including 
the maxima at fixed m^/2) significantly differ. The largest differences are found at low 
mi/2 and at low tan/3. 

A variety of numerical issues may cloud the evaluation of the ^ parameter in super- 
symmetric models connecting the GUT scale to the weak scale. For example: 

• One problem is that the convergence at low values of fi will depend in part on the 
initial guess for the supersymmetric masses at the very beginning of the iterative RG 
evolution process. The default guess for the supersymmetric masses used in the two 
codes respectively reads 

msusY = Y^mfT^M^ Isajet 7.72 (2.8) 
msvsY = 0.5 (mo + M1/2) + mz SUSPECT 2.34 (2.9) 

Clearly, the different trial values used by the two codes can alter the final convergent 
solution and the value of fi as well, particularly in the highly fine-tuned region where 

• A further problem occurs in that SUSPECT 2.34 uses two-loop RGEs only for gauge 
and Yukawa couplings, but not for soft SUSY breaking terms, while Isajet uses two- 
loop RGEs throughout its RG treatment. 

• Another difficulty occurs in that loop corrections must be included in the formulae 
for the minimization of the scalar potential of the theory. These loop corrections 
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Figure 4: A comparison, at mt = 174.3 GeV and along iso-mi/2 slices, of the ^ parameter as a 
function of mo between the two RG evolution numerical codes Isajet 7.72 (red dashed lines) 
and SUSPECT 2.34 pi (black solid lines), at tan/3 = 10 (left) and 50 (right). 



depend on the full spectrum of supersymmetric particles. However, to calculate the 
full spectrum, the value of fi must be known. It is possible in the low /u region that 
the tree level value of fi"^ < 0, while loop corrections will lift /x'^ > 0. In this case, a 
guess must be made as to what the loop corrected value of n is, just so that a viable 
spectrum of SUSY particles can be calculated, and used as input, in its turn, to the 
loop corrections. The value of fi will depend on how this guess is made. 

• As evident from Eq. (^]^), at very large values of the common scalar mass, is 
extremely sensitive to the top Yukawa coupling ft. The latter is defined as 

ft = V2 mt/vu, (2.10) 

where mt is the running top quark mass in the DR scheme. The computation of rrit 
suffers from uncertainties related to (1) the extraction of the DR top mass from its 
pole or MS values (i.e. the inclusion of Standard Model threshold corrections) and 
(2) the implementation of the SUSY loop corrections, and the scale choice at which 
these are implemented. Similar ambiguities pertain to the evaluation of the bottom 
and tau Yukawa couplings. 

• The numerical results will depend on the scale choice at which different SUSY parti- 
cles are integrated out of the effective theory. A related issue is the choice of scale at 
which various SUSY threshold corrections are implemented. For instance, SoftSUSY 
and Spheno assume the MSSM is valid from Mqut all the way to Mz- When scalars 
have masses of several TeV, such as in the HB/FP region, this may not be a good 
assumption. 

As a bottom line, large numerical uncertainties plague the study of the phenomenology 
of the HB/FP region at the GUT scale, depending on a number of assumptions in the 
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details of the RG evolution which lead to significant discrepancies in the determination 
of /i from REWSB; on top of this computational ambiguity, a possibly larger uncertainty 
stems from the input value of the top quark mass: even at the level of accuracy with which 
rrit will be measured at the LHC the induced variations in the low-scale parameters in the 
HB/FP region would give rise to completely different phenomenological scenarios at the 
same GUT-scale input values [20|. Lastly, currently available numerical codes do not fully 
access the phenomenologically interesting mSUGRA region at ^ <C Mi, which can thus be 
only explored resorting to a low-scale parameterization. 



2.2 Outline of the low-energy parameterization 

We described above how the procedure of outlining a systematic mapping between the phe- 
nomenologically relevant low-energy parameter space and the customary mSUGRA GUT- 
scale input setting faces a number of intrinsic issues. Moreover, those issues evidently 
blur the possibility of a bottom-up reconstruction of mSUGRA high-energy input param- 
eters. Our main point here is therefore that the collider and dark matter phenomenology 
of the focus point region is better studied through a convenient two-parameters low-energy 
description, which is complementary to the usual GUT-scale parameterization, and which 
thoroughly reproduces all phenomenological implications of the "mother theory" at the high- 
energy scale. This section is henceforth devoted to the construction of such a low-scale 
parametrization, of which we will then make use in the remainder of this report. 

We argued above that the physical parameter space of the HB/FP region is given 
by the set (Mi, //, tan/3). At a given value of tan/3. Mi (and therefore of mi/2 at the 
GUT scale) and of the input top quark mass rut, REWSB bounds the range of /i from 
above, through the condition given in Eq. ( |2.4D . As pointed out in the previous section, 
the precise boundaries of the region allowed by REWSB will also depend on numerical 
subtleties. However, we will show below that the REWSB boundary at large ^ » Mi 
lies outside the cosmologically relevant HB/FP region. As a result, these uncertainties 
will not affect the discussion of the HB/FP region phenomenology, since that portion of 
the parameter space will be disregarded after the neutralino relic density analysis. The 
other soft breaking gaugino masses M2 and M3 are given, as functions of Mi, by the GUT 
relations specified in Eqs. (p.5D -(^). 

Since the heavy scalar sector is largely decoupled from the low-energy phenomenology 
in the HB/FP region, the details of the sfermions and heavy Higgses spectrum play a very 
marginal role. In this respect, we resort in most of our plots, to setting, for simplicity, all 
sfermions and heavy Higgs masses to a common value mj = 5 TeV, much larger than all 
other relevant weak-scale parameters, and close to the average value of the masses obtained 
by a full RGE treatment. Further, we set all trilinear couplings to zero. 

The value of the lightest Higgs mass enters, instead, quite critically in a few 
quantities, in particular those related to the neutralino scattering off matter: for instance, 
the neutralino-quark spin independent scattering cross section uf^ , in the limit of large 

Ziq 

scalar masses, is dominated by i-channel light Higgs exchanges, and scales as crf^ oc 

Ziq 

m7^. In mSUGRA, the value of the lightest Higgs mass depends, in principle, on all 
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input SUSY parameters. In the HB/FP region, at large universal scalar masses, we find 
a critical dependence on 777-1/2 > while the GUT-scale values of mg and ^0 are far less 
important. Consistently with the assumption of a common low-scale scalar mass mj, we 
therefore expressed m/^ as a function of the parameters Mi and tan/3, according to a 
phenomenological functional dependence of the form 

mH = ml (^y. (2.11) 

The formula was then fitted against mSUGRA points randomly generated at large values 
of rriQ. We find that and m^j do not critically depend on tan (3, and the best fit values 
read = 115.5 GeV and = 100 GeV, respectively. The exponent 7 depends instead 
more sensitively on tan/3: for instance, we find 7 ~ 0.03 for tan/? = 10, and 7 ~ 0.0335 
for tan/3 = 50. We estimate the typical accuracy of Eq. ( 2.11| ) in reproducing the correct 



value of rrih to be less than 1 GeV for Mi > 300 GeV, and within 2 GeV for values of Mi 
smaller than 300 GeV. The relative error induced in cr~^ is therefore expected to be of 
0(1%). 

2.3 Neutralino relic abundance 

Large values of the sfermion masses in the HB/FP region help alleviate a number of well- 
known phenomenological difficulties of supersymmetric extensions of the standard model, 
ranging from the SUSY flavor and CP problems, to dimension-five proton decay operators 
which often appear in SUSY-GUT embeddings. In the present framework, where we neglect 
CP violating phases, and assume a minimal flavor structure, the SUSY contributions to 
rare processes, e.g. the branching ratios b 57, Bg — > fj,^fi~, or to the muon anomalous 
magnetic moment, are very suppressed, being mediated by sfermions loops, and do not 
constrain the HB/FP region parameter space. Further, limits from direct sfermion searches 
at colliders are always very distant from the sfermion mass scales of the HB/FP region. 

The phenomenological constraints which apply to the HB /FP region of mSUGRA are 
henceforth limited to the LEP2 searches for the lightest chargino, and to gluino searches 
at the Tevatron. In the present setting with universal gaugino masses, the latter bound is 
however much less constraining than that stemming from chargino searches. We use here 



the mass limit m^y^ > 103.5 GeV ||2^, although in the pure-higgsino region {fi <C Mi) the 
LEP2 bound is weakened by the quasi-degeneracy between the lightest chargino and the 
LSP (namely, n^'^/^ > 92.4 GeV |2^). The (Afi,/x) plane is further constrained, at a given 
value of tan /?, by the REWSB conditions which limit the range of fi from above. In the 
left-over portion of parameter space we compute the thermal LSP relic abundance O h? 
with the DarkSUSY package [29|. We rule out models giving ^^i^"^ ^ 0.13, according to the 



2-cr upper bound on the cold dark matter abundance derived by the WMAP collaboration 
In Fig. 1^, left, we plot the parameter space of the HB/FP region using the scheme 



m 



outlined in Sec. 2^, at tan/? = 50 and positive fj,. In this plot only, the scalar masses are 
determined by the value of mo needed to give the appropriate jj, value as in Fig. ^6. The 
region shaded in red is excluded by LEP2 chargino searches, while in the gray shaded area 
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m,„ (GeV) 




Figure 5: (Left): The allowed region, on the (Afi,/i) plane, at tan/3 = 50. The area shaded in 
red features a lightest chargino mass below the LEP2 limit |Q, i.e. < 103.5 GeV. In the gray 

shaded area the value of fi exceeds the maximal one allowed by radiative electroweak symmetry 
breaking (see Fig. 0). The yellow shaded area is cosmologically excluded since the thermally pro- 
duced neutralino relic density exceeds the 2-a upper bound on the cold dark matter abundance. 
Finally, points in the green region feature a neutralino thermal relic abundance compatible with 
the upper and lower limits on the cold dark matter content of the Universe. We also indicate the 
location of the coannihilation strip in the (Afi, /i) plane. The solid black lines indicate those points 
giving exactly the central value ^^i^^ ^ ^-^^ HB/FP region. (Right): The neutralino ther- 

mal relic abundance as a function of /i at fixed Mi = 250 GeV (solid black line), 1000 GeV (red 
dashed line) and 1500 GeV (blue dot-dashed line). The green strip indicates the WMAP favored 
range for the cold dark matter abundance. 



the value of fi exceeds the maximal value compatible with REWSB at tan /? = 50 and 
/i > (see Fig. |l|, right). The parameter space portion shaded in yellow indicates where 
the upper bound on ^^_^h'^ is violated. The viable parameter space of the theory is thus 
restricted to the green band (giving r^^^/i^ within the 2-a WMAP range) and to the region 
in white {Q^^h^ below the 2-a WMAP range). In the low-relic density white region we 
suppose that some mechanism in the Early Universe has enhanced the final relic density of 
neutralinos {e.g. through non-thermal neutralino production |31], or through a modified 
cosmic expansion at the time of neutralino freeze-out, as it might be the case in scenarios 
with quintessence ||3^, with a Brans-Dicke- Jordan modified theory of gravity |32, 34 1 or 
with an anisotropic primordial expansion of the Universe p3| , p5|]). We therefore work 
under the assumption that all CDM is composed of relic neutralinos. 

For completeness, we also indicate in Fig. |5| the location of the stau coannihilation strip 
in this low-energy representation, squeezed along the line of maximal ^ values. Had we 
chosen to consider negative values for /i, a rapid heavy Higgs s-channel exchange-mediated 
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Figure 6: The iso- level curves for the lightest neutralino mass (left) and composition (in terms 
of the higgsino fraction = -/V13 + -/V^4, right), on the parameter space outlined in Fig. |[ at 
tan/3 = 50 and > 0. 



annihilation funnel would also have opened. The latter would lie at rather large values of [i 
(slightly below the coannihilation strip, since the corresponding mg would have been larger 
in the funnel than in the coannihilation region, and therefore the corresponding range 
would be slightly lower (see Fig. possibly extending to even larger values of M\. In this 
respect, the (Mi,/x) plane is not a good representation of the coannihilation and funnel 
regions of mSUGRA, since it trades a physically relevant parameter for those regions (?7io) 
for a parameter which is instead not critical for il^^/i^, i.e. fi. We therefore stress a strong 
complementary role of the GUT-scale and of the present low-scale parameterizations for 
the phenomenological study of the mSUGRA model. 

The shape of the parameter space region giving exactly the central value of the CDM 
density determined by WMAP is indicated with a black line approximately lying in the 
center of the green-shaded 2-a area. The shape of that line provides a non-trivial infor- 
mation on the mechanisms responsible for a WMAP-compatible neutralino thermal relic 
abundance in the HB/FP region. In the large mi/2 (large Mi) limit, the neutralino relic 
abundance corresponds to that of a pure higgsino-like LSP, and is fitted, with high preci- 
sion, by the formula 

/ u \ 1.906 

/i2~0.l(— ^) . (2.12) 

^1 VlTeV/ ^ ^ 

This formula entails the important result that the maximal neutralino mass compatible at 
95% C.L. with the WMAP upper limit on the CDM abundance in the minimal supergravity 
model corresponds to 1150 GeV (We recall that in the funnel and coannihilation regions 
is always found to be less than 900 GeV, see e.g. [^]). On the other hand, the central 

WMAP CDM abundance value here is attained at m;? ~ 1050 GeV. 

^1 

The central part of the curve lies along the Mi ~ /i edge; in this region the interplay of 
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coannihilation processes and of a mixed bino-higgsino LSP cooperates to fulfill the condition 
Q^^/i^ ~ ficDM^^- In the (Mi,^) plane the 2-(j range of n shrinks: this fact depends on 
the transition from the pure higgsino to the pure bino regime. We illustrate explicitly this 
phenomenon in the right panel of fig. |5|, where we show ^^Ji^ as a function of ^ along 
slices at fixed Mi. The mentioned transition, which we indicate as "coannihilation regime" 
in the figure, starts at ^ ~ Mi and ends when ^ ^> Mi; in this regime, the relic abundance 
is clearly a very steep function of //. 

As Ml is further decreased, the neutralino mass compatible with the WMAP CDM 
relic abundance decreases, while the bino fraction and the mass splitting between the LSP 
and the lightest chargino increase. This is illustrated in Fig. which also shows the iso- 
neutralino mass contours and the iso-higgsino fraction levels on the {Mi,fi) plane, which 
will be the object of our phenomenological analysis in the remainder of this report. 

In particular, the right panel of Fig. ^ illustrates that for mi/2 ^ 2.5 TeV the lightest 
neutralino in the HB/FP region is always a very mixed bino-higgsino particle, with higgsino 
fractions Zh = N^^ + N^^ between 0.1 and 0.9. At larger mi/2, instead, Zh > 0.9, and the 
"pure higgsino regime" {Zh > 0.99) is reached at mi/2 ^ 3 TeV. 

Finally, let us point out the threshold effects, at Mi ~ 200 GeV and Mi ~ 100 GeV, 
respectively, corresponding to m^^ ~ mt and m^^ ~ mw,z- When one of the tt, WW, ZZ 
final channels closes up, the resulting suppression of the lightest neutralino pair annihilation 
cross section must be compensated with an increased higgsino fraction and a reduced mass 
splitting between the LSP and its coannihilation partners, resulting in the above mentioned 
bumps at the corresponding thresholds. 

As a concluding remark, we stress that the picture we outlined above is very mildly 
dependent on the particular value of tan P we picked, and on the sign of /U as well. We 
explicitly worked out the WMAP allowed parameter space for lower tan P, and the emerging 
picture is almost indistinguishable from what we show here, providing evidence for a kind 
of universality in the HB/FP region phenomenology on the (Mi,//) plane. In particular, 
our conclusions on the maximal LSP mass in mSUGRA are not significantly affected. 



3. Dark Matter Phenomenology 



In the present section we study the implications of SUSY models belonging to the HB/FP 



region parameter space outlined in Sec. 2.3 for dark matter searches. Following previous 
analysis |37, we make use of two extreme halo profiles, a cored halo model (the 

Burkert profile, see ||4^, 41, ^) and an Adiabatically contracted version of the cuspy 
halo model proposed in Ref. ||43] (which we dub Adiabatically Contracted N03 profile). We 



claim that those two instances are indicative of the range of variations in dark matter 
detection rates induced by different consistent models of the dark matter distribution in 
the halo, respectively, giving minimal (Burkert) and maximal (Adiab.Contr.N03) rates. 
We refer the reader to Ref. [pTl for details. 



Sec. 3J. is devoted to a parallel analysis, for the two halo models, of the current con- 
straints coming from antimatter and gamma ray fiuxes. We also impose the constraint 



coming from the overproduction of Li in the Early Universe |45, 4(:]. Sec. 3.2 is instead 
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Figure 7: Current constraints on the (Mi,/j,) plane from dark matter physics, at tan/3 = 50, 
/i > 0, for the Burkert Halo Model (left) and for the Adiabatically contracted N-03 Halo Model 
(right). The gray shaded area is excluded by the LEP2 searches for the lightest chargino and for an 
excessive neutralino thermal relic abundance (see Fig. ^), while the green shaded region is the 2-cr 
WMAP-favored region from the neutralino relic abundance. The area shaded in dark blue is ruled 
out by current data on the antiprotons flux, while the area shaded in brown by the primordial ^Li 
abundance bound. In the plot to the right, the cyan region is not consistent with positrons flux 
data, and the region shaded in black gives an excessive gamma-ray flux in the energy range of the 
EGRET experiment. 



devoted to prospects for dark matter detection in the HB/FP region at future experi- 



ments, while in Sec. |3.3| , we study the spectral features at future space-borne antimatter 
experiments of models in the HB/FP region giving a thermal relic density of neutralinos 
consistent with the WMAP CDM abundance. 



As mentioned in Sec. 2^, we assume here that models with a low thermal rehc abun- 
dance, in the /i <C Mi region, are responsible for all the CDM in the Universe, in virtue of 
the above mentioned mechanisms of relic density enhancement (see e.g. Ref. ||35[| ). 

3.1 Overview of current constraints 

Current dark matter detection results include the recently released direct detection exclu- 
sion limits delivered by the CDMS collaboration ||4^, the SuperKamiokande upper limit on 
the muon flux from neutralino pair annihilation in the core of the Sun ||4^ , the antimatter 
fluxes as measured by balloon-borne experiments and the EGRET data |5C] on the 



gamma ray flux from the Galactic Center. Besides direct and indirect dark matter searches, 
neutralinos are also constrained by the synthesis of ®Li in post-freeze-out neutralino anni- 



hilations, as recently shown in Ref. [45|. This constraint can be rephrased as a constraint 



on the neutralino pair-annihilation cross section |45], and, contrary to dark matter search 



results, is independent of the halo model under consideration. 
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In Fig. 0, we depict the above mentioned current exclusion limits in the (Mi, ^) plane. 
The region shaded in gray corresponds to the parameter space excluded either by LEP2 
searches for charginos, or by the lack of REWSB solutions, or by the overproduction of 
neutralinos in the Early Universe. The green-shaded area corresponds to the models giving 
^Zi^'^ ~ f^CDM^'^ at tlie 2-0" level. The left panel corresponds to the Burkert profile, while 
the right panel to the Adiab.Contr.N03 profile. 

The ^Li constraint is violated on the brown-shaded area, extending in the pure higgsino 
limit in the mass range 90 GeV < ^150 GeV, where the neutralino pair annihilation 
cross section {crv)Q is maximal. 

Neutralino induced primary antimatter fiuxes, derived according to the approach de- 
scribed in detail in Ref. |^^, are not consistent, at 95% C.L., with the measured total 
(signal plus background) antimatter flux of antiprotons and positrons [^] in the regions 
shaded in dark and light blue, respectively. We notice that for both halo models, the 
antiproton fiux constraint is stronger than the ^Li bound; the positron fiux is instead in 
excess to the measured one only for the cuspy halo model, at m^^ < 170 GeV. 

The EGRET data give an energy-dependent upper bound on the gamma ray flux from 
the galactic center, under the conservative hypothesis that the neutralino pair-annihilation 
induced signal dominates over a negligible background. The gamma ray flux depends 
critically on the inner structure of the halo profile, and is therefore greatly enhanced for 
a cuspy profile, while it is suppressed for a cored halo model. In fact, the EGRET data 
do not give any constraint if one assumes the Burkert profile, while they rule out pure 
higgsinos as heavy as 350 GeV with the cuspy Adiab.Contr.N03 profile. 

As a final remark, we point out that neither direct detection searches nor current data 
on the muon flux from the Sun give any constraint on the HB/FP region parameter space. 



3.2 Future dark matter searches 

We outline in Fig. |^ (for the Burkert profile) and ^ (for the Adiab.Contr.N03 profile) the 
sensitivity reach of a number of future direct and indirect dark matter detection experi- 
ments. Models in the area enclosed within each contour give a signal which is larger than 
the expected sensitivity of the corresponding search facility. 

"Stage-2" detectors refer to experiments like CDMS2 Q, Edelweiss2 [||], CRESST2 H, 



ZEPLIN2 1 54], which will be operative in the near future (the reference sensitivity curve we 
take here is that corresponding to the CDMS2 experiment pTI). We indicate as "Stage-3" 



detectors ton-size experiments like XENON [|5|, Genius pi, ZEPLIN4 |5l and WARP p| 



(the reference sensitivity curve is here chosen to be that of XENON [|55[ ) . 

The reach of Stage-2 detectors is found to be limited to a tiny region at very low 
masses, largely already ruled out by the chargino mass limit. A neutralino in the HB/FP 
region is therefore not likely to be detected at Stage 2 direct detection experiments. On 
the other hand, ton-sized detectors look very promising, their sensitivity extending, quite 
independently of the halo model, and largely independently of tan /3, over the whole re- 
gion compatible at 2-a with the WMAP CDM neutralino thermal relic abundance, up to 
^1/2 ^ 3 TeV. At larger values of the common gaugino mass, the higgsino fraction becomes 
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Figure 8: Reach contours of future dark matter search facilities. The region shaded in gray is 
excluded in the light of the results shown in Fig. || and ^ The value of tan/3 is set to 50, /i > 0, 
and the cored Burkert Halo Model is assumed. 



exceedingly large (see Fig. ^, right), and, since o"|^p oc N^{1 — N^)'^, the resulting direct 



.g...;, ...u, ...... u^' - ^r2n 

detection rates are suppressed. 

The flux of muons from the Sun is also extremely sensitive to the degree of gaugino- 
higgsino mixing. The equilibrium between the capture rate inside the Sun and the neu- 
tralino pair annihilation, basic in order to produce a sizable neutrino flux out of the Sun, 
is reached provided the neutralino features (1) a sufficiently large pair annihilation cross 
section and (2) a large enough spin-dependent neutralino-proton scattering cross section. 
These two conditions explain the shape of the IceCube reach contours [^], which extend 
along the largest {crv)o area at low fi as far as the gaugino fraction is still non- negligible, and 
in the maximal mixing region. The latter region largely overlaps the WMAP favored green- 
shaded region. Remarkably enough, neutralinos producing the WMAP required amount of 
relics to make up the CDM in the Universe tend to maximize, in the HB/FP region, the 
detection rate at neutrino telescopes! We find that, quite independently of the halo model, 
the IceCube reach along the cosmologically favored strip extends up to mi/2 ^ 1-5 TeV, 
corresponding to a neutralino mass between 600 and 700 GeV. 



The dependence of the antimatter flux on the halo model |3£] is found, instead, to be 
indeed critical. For antiprotons and positrons, the largest fluxes correspond to the pure 
higgsino region, partially extending into the large bino-higgsino mixing region. The shape 
of the reach contour for the Pamela experiment (computed following the approach outlined 
in Ref. ||3^) reflects the role of the top threshold (below which a much larger higgsino 
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Figure 9: Same as in Fig. but for the Adiabatically contracted N-03 Halo Model 



fraction is needed) and extends into the light neutralino mass region until the second 
critical threshold {ttl^^ < mw) is reached. Assuming a cored profile, models giving the 
WMAP relic neutralino abundance do not yield large enough antiprotons and/or positrons 
fluxes, while with a cuspy profile the reach along the WMAP strip in the antiproton channel 
extends up to neutralino masses as large as 400 GeV. 

Low energy antideuterons have been shown to provide a clean indication of new physics, 
since the background flux in the 0.1- 1 GeV antideuteron kinetic energy interval is ex- 
tremely suppressed pO|. We assessed the antideuteron flux for the AMS-02 experiment 



|61], which will be sensitive to a flux in the energy band 0.1 — 2.7 GeV at the level of 
4.8 X 10^^ m^^ sr^^ GeV~^ sec~^, and for a GAPS experiment [^2| placed on a high lati- 
tude mission satellite, sensitive to antideuterons in the energy band 0.1 — 0.4 GeV at the 
level of 2.6 x 10^^ m^^ sr^-*^ GeV~"^ sec~^. We find that antideuteron fluxes accessible to 
AMS-02 are in general excluded by current bounds on the antiprotons flux, independently of 
the assumed halo profile. The GAPS sensitivity extends instead well beyond the antiproton 
search reach of Pamela in the WMAP favored region, with maximal accessible neutralino 
masses as large as 400 GeV in the case of a cored profile and of even 700 GeV in the case 
of a cuspy profile. The antideuteron flux is also greatly sensitive to the gaugino-higgsino 
mixing, thus, as in the case of the muon fiux from the Sun, the maximal reach is gained 
exactly along the WMAP 2-cj region. 

The parameter space reach of the GLAST experiment ||6^ , finally, largely depends on 
the details of the inner structure of the halo model under consideration. As can be deduced 
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Figure 10: The signal-over-background for antiprotons (left) and positrons (right), as functions of 
the antimatter particles' kinetic energy, for four values of the neutralino mass along the WMAP- 
relic-abundance line at tan/3 — 50, /i > 0. See the text for details on the background and signal 
computation. 



comparing Fig. |8| and ^, the resulting GLAST reach can be either the worst or the best 
among all direct and indirect detection channels, depending on the amount of dark matter 
in the center of the Galaxy: for the two halo models under consideration, for instance, this 
induces a variation of more than four orders of magnitude [^]! Little can therefore be said 
about the sensitivity reach of future gamma rays experiments without strong astrophysical 



assumptions (for recent related studies see 1 64 , 14 , ) 



3.3 Spectral features at space-borne antimatter searches 

In the preceding section we assessed the reach of future space-based antimatter experi- 
ments by means of a statistical treatment which evaluates the possibility of disentangling a 
pure-secondary antimatter flux from the presence of a statistically non-negligible primary 
component. In case such a signal is detected at Pamela or AMS-02, we will be given the 
opportunity of studying in some detail the spectral features of the primary component, 
depending on the relative signal-to-background {S/B). In particular, the future wealth of 
data on antimatter fluxes will greatly reduce the uncertainties in the background determi- 
nation to an unprecedented level of accuracy. For this reason, we investigate in this section 
the S/B yields of models in the HB/FP region, concentrating on the parameter space slice 
giving ^^_^h'^ equal to the central value of the WMAP-inferred CDM abundance Ocdm^^- 



We show in Fig. 10 the S/B as a function of the antiparticle's kinetic energy for 
antiprotons (left) and positrons (right), at four different neutralino masses. The antiprotons 
S/B features in all models a clean peak at Tp ~ O.lm^^, generated by hard decay modes, 
particularly from gauge boson decays, plus a low energy tail, mainly fueled by products of 
bb final state processes. Positrons show a more complex S/B, featuring a series of peaks. 
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Figure 11: On the right y-axis of the two panels we plot the antiparticles' kinetic energy Tmax at 
which the maximal {S/B) is attained, in the case of antiprotons (left) and positrons (right). The 
left y-axis indicate the actual maximal {S/B) at the corresponding T,„ax- In the case of positrons 
(right panel), we also plot the {S/B) at = m-^^/2. 



corresponding to quark jets yielding vr"*" — > e"*" and positrons from r"*" decays at low energies, 
and to positrons generated by gauge boson decays at higher energies. In particular, prompt 
W and Z decays motivate the bump at ~ m^^/2, which is however very suppressed, 
and hard to recognize, particularly beyond the = rrit threshold, which tends to soften 
the positron spectrum. 

The location of the maximal {S/B) ratios is studied in Fig. 11, where we plot on the 
right axis the corresponding antiparticle's kinetic energy, and on the left axis the actual 
S/B value at the maximum. In both cases, the mt threshold (corresponding to Mi ~ 200 
GeV in the plot) is clearly visible. The maxima, in that low neutralino mass end, are 
located around 10 GeV, for both positrons and antiprotons. For m^^ > mt, the location 
of the antiproton's maximal S/B approximately linearly tracks m^^, while for positrons 
the maxima are positioned around 10 GeV, until eventually the S/B corresponding to the 
Te+ ~ m^^/2 prompt positron production dominates. 

In order to understand to which extent the above analyzed spectral features are spe- 
cific to the HB/FP region (or to any mixed higgsino-bino LSP scenario), we carried out a 
comparison of {S/B) in different LSP scenarios, namely that of a pure bino (mainly anni- 
hilating into bb and r+r^ pairs), as in the coannihilation and funnel region of mSUGRA, 
and that of a pure wino LSP (mainly annihilating into W^W~ pairs), as in the minimal 
anomaly mediated supersymmetric breaking model. The case of the bino is characterized 
by a much larger antimatter yield at low energies, which, for instance, smooths out the 
maximum in the antiproton S/B pointed out above, and which shifts the positron's maxi- 
mal S/B towards much lower energies. This LSP scenario is thus virtually distinguishable 
from the HB/FP region LSP scenario on the basis of a correlated S/B analysis. The case 
of a wino-like LSP is instead more subtle. The antiproton S/B has a more depressed 
low-energy tail (which can be however hard to disentangle due to uncertainties in the low- 
energy antiproton background computation), and the location of the maximal S/B as a 



-19- 




M| (GeV) 

Figure 12: mSUGRA with tan/3 = 50, /i > 0. (Left): The mass spectrum in the neutrahno- 
chargino sector, along the same parameter space line. (Right): The higgsino fraction, as a function 
of Ml, along the line at h"^ = 0.11 (red dashed line), and the absolute mass splitting between 
the lightest chargino and the lightest neutralino along the line at ^^J^"^ = 0.11 (black solid line). 



function of the mass is at larger Tp/m^^ values. On the other hand, a possible handle is 
provided by a much stronger peak in the positron spectrum corresponding to Tg+ ~ m^j2. 

The occurrence of a maximum in the antiproton's {S/B) correlated with the neutralino 
mass can evidently be used as an indirect indication of the neutralino mass scale, keeping 
in mind the above mentioned caveat on a possible entanglement of the mixed-higgsino- 
bino and wino LSP scenarios. Vice versa, should collider experiments or other dark matter 
searches point at a certain neutralino mass, the analysis of the antiproton's S/B can be 
used as a cross-check to understand the nature of the dark matter particle. Positron fluxes 
appear, instead, less promising for the same task, the absolute values of the S/B being 
moreover much smaller than those of antiprotons. The exciting perspective of a correlation 
between the antiproton's and the positron's spectral features, which could point to a specific 
LSP dark matter scenario, therefore also appear rather problematic. 



4. HB/FP region at LHC and ILC 
4.1 HB/FP region at the CERN LHC 

The reach of the LHC in mSUGRA's (mo, ^1/2) plane was described in Ref. |]l8| (see |17| 



for related earlier work). The search strategy was based on the detection of gluino and 
squark cascade decay products, namely multiple jets and/or leptons and large transverse 
missing energy. Since sfermion masses in the HB/FP region lie in the multi-TeV range, 
LHC will not produce squarks and sleptons at detectable levels. Gluinos can be relatively 
light {rrig < 2 TeV, left frame of Fig. |l^) only in the low Mi ~ 250 GeV (or, equivalently, 
low ij) part of the HB/FP region along the line ^^_h'^ = 0.11. It was found in [18| that 
in the HB/FP region of mSUGRA, gluino masses of up to ~ 1.8 TeV could be probed. 
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Figure 13: The Wi — Z\ (left) and — Z\ (right) mass spUttings, as functions of A/i, for four 
parameter space shces, at tan/3 = 5, 50 and > and < 0, along the lines giving ^^_^h^ = 
^CBuh^- We also include gray points obtained via a random scan over the HB/FP region of 
parameter space. 



Therefore one would not expect the method, based on the production of heavy strongly 
interacting superpartners, to work very well in the HB/FP region for larger Afi values. We 
note that recent work by Mercadante et al.^^ employed 6-tagging in their analysis and 
have achieved an up to 20% extension of the LHC reach in the HB/FP region. 



However, as can be seen from the left frame of Fig. 12, charginos and neutralino 
(collectively dubbed as -inos) can still be relatively light in the HB/FP region. Even more 
encouraging is the fact that the mass splitting between the lightest chargino and lightest 
neutralino stays > 20 GeV all the way up to Mi ~ 850 GeV and the mass gap never exceeds 
the Z boson mass (right frame of Fig. |T^ ). This means that 3-body decays Wi — > Ziff and 
Z2 Ziff are open, and one expects multiple leptons in the final state. We will resort to 
the leptonic signals when exploring -ino production at the LHC, since soft jet signals have 
an enormous QCD background. 

The Wi — Z\ {Z2 — Zi) mass gap is shown in the left (right) frame of Fig. |l^ for various 
values of tan/3 and signs of ^. Positive fi is beneficial for our purposes, because the mass 
gap is larger than for negative /U, and therefore the leptons in the final state are harder. 
A positive value of /i is currently favored by the discrepancy between the measured muon 
g — 2 value and the one calculated in the SM using e^e~ data for the diagrams involving 



hadronic vacuum polarization [36|. 



For the rest of this section we restrict ourselves to the following values of input param- 
eters: tan/3 = 50, ^ > and mt = 174.3 GeV. Unless stated otherwise, we will be working 
along the WMAP favored ^^j^"^ ~ ^-^^ ™ HB/FP region. Along this line Mi and 
fj, are not independent, and one needs to specify only the Mi value. 
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Figure 14: In the upper frame, we show the total cross section and the cross section after the 
pre-cuts {EJp^'^ > 200 GeV and at least 2 jets with pr > 40 GeV and \-q\ < 3) for SUSY particle 
production at the LHC along the VlcDMh^ = 0.11 line, as a function of Mi. In the lower frame, 
we show LHC cross sections after optimized cuts for various gluino cascade decay event topologies 
and corresponding 5cr discovery limits for 100 fb^^ integrated luminosity along the WMAP favored 
region. LHC can probe to Mi ~ 250 GeV, which corresponds to mg ^ 1.8 TeV. 



4.1.1 Reach via gluino cascade decays 



We use Isajet 7.72 [23| for the simulation of signal and some of the background events at the 
LHC. A toy detector simulation is employed with calorimeter cell size A?] x = 0.05 x 0.05 
and —5 < ?7 < 5. The hadronic energy resolution is taken to be 80%/\/^ for |r/| < 2.6 
and 100%/\/^ for |??| > 2.6. The electromagnetic energy resolution is assumed to be 
3%/^/E. We use a UAl-like jet finding algorithm with jet cone size R = 0.5 and p^r^^ > 25 
GeV. Leptons are considered isolated if the visible activity within the cone AR < 0.3 is 
Y^^^lis ^ 2 GeV. The strict isolation criterion helps reduce multi-lepton background from 
heavy quark (especially tt) production. Leptons (es or /is) have to satisfy the requirement 

require that leptons would have \r]\ < 2.5 and jets would 

propagate within |?7| < 3. 



First, we re-plot the reach of the LHC in Fig. 14, using the procedure described in |18 



All events had to pass the pre-cuts, which impose the requirement that i?™*'^* > 200 GeV 
and there are at least 2 jets with p^^* > 40 GeV. The definitions of jets and leptons can 
be found in |18], as well as the description of the cut optimization procedure. We choose 
the cuts which were found to be optimal for the HB/FP region. The events are divided 
into several classes, characterized by the number of leptons in the final state. The upper 
frame of Fig. 14 shows the total supersymmetric particle production cross section and the 
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Figure 15: The leading order production cross sections at the LHC for all possible pairs of neu- 
tralinos (left), and of a chargino plus one of the two lightest neutralinos (right), along the line at 
= 0.11. 

^1 



cross section after pre-cuts. The lower frame of Fig. 14 shows the signal cross sections after 
the optimal cuts in various channels. The 5a discovery reach for 100 fb~^ of integrated 
luminosity is shown by short horizontal lines for each channel. One can see that the cuts, 
which are optimized for gluino pair production and subsequent cascade decays, provide a 
reach by the LHC of up to Mi ~ 250 GeV, corresponding to a value of rrig ~ 1.8 TeV. At 
that point the total sparticle production cross section is still sizable: for Mi = 250 GeV, 
the total SUSY cross section is ~ 400 fb. This fact motivates us to look next at -ino pair 
production at the LHC. 

4.1.2 Trilepton production at the LHC in the HB/FP region 

Total production cross sections for neutralino pair production at the LHC along the line of 



constant ^^Jt^ = 0.11 are shown versus Mi in the left frame of Fig. 15. The right frame 



of Fig. 15 shows total cross sections for associated production of the lightest or next-to- 



lightest neutralino with a chargino. Similarly, we present the associated production of 3rd 



or 4th lightest neutralino with a chargino in the left frame of Fig. |16| and the chargino 
pair production cross section in the right frame of Fig. |l6|. Of all the -ino production cross 
sections, WiZi, W1Z2 and W^W^ production are generally the largest. 

The gluino pair production cross section is presented in Fig. ^ versus Mi . We have 
also shown the total -ino production cross section for comparison. Assuming 100 fb~^ 
integrated luminosity, LHC would produce less than 10 gluino pairs for Mi ~ 300 GeV, 
while ~ 10^ -ino pairs would be produced. 

Given the relative total production cross section rates, it might be beneficial to exam- 
ine -ino pair production signals in the HB/FP region, as well as gluino pair signals. Single 
lepton signals from W\Z\ where W\ — > lviZ\ will be buried under an immense background 
from direct W boson production. Likewise, dilepton production from reactions such as 
W^j^W-i" production will be buried beneath large backgrounds from W~^W~ and produc- 
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Figure 16: The leading order production cross sections at the LHC for aU possible pairs of a 
chargino plus one of the two heaviest neutralinos (left), and of two charginos (right), along the line 

at = 0.11. 

^1 



tion. Four lepton signals from reactions such as Z2Z2 production would be very distinctive 
and it is possible to find the cuts which reduce the SM background. However, our prelim- 
inary analysis found that the 4 lepton signal rates fall very quickly with increasing Mi, 
and this channel did not provide any additional reach compared to the optimized cuts for 
gluino pair production. Thus we were led to consider the clean trilepton signature at the 
LHC in more detail [|67|]. This signal provides the best reach in mSUGRA at the Tevatron 
1 16, due to W1Z2 production and subsequent decays Z2 — > II Zi and Wi — > IviZi. 

Let us first examine the dominant production processes, which could produce 3 or 
more leptons in the final state, at Mi = 250 GeV, where the reach due to optimized cuts 
peters out: 



a{WiZ2) = 70.8fb, 
a{WiZs) = 76.8/6, 
a{Z2Z3) = 28.0/6, 
a{W2W2) = 17.2/6, 
(7(W2^4) = 36.4/6. 
The relevant branching fractions are: 
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where / stands for either e or fi. It is now possible to estimate the maximal possible 3 
lepton event rate before any cuts for 100 fb~^ integrated luminosity at the LHC. We do 
that only for the processes (pi]), and (|43| ) here: 

N3iiWiZ2) = 109, 
NsiiWiZs) = 118, 
N3liZ2Z3) = 14. 
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Next we proceed to the fast sim- 
ulation. The dominant backgrounds 
for the clean trilepton signature are tt, 
WZ and ZZ production. When eval- 
uating WZ production, it has been 
shown to be of crucial importance to 
evaluate the full qq' iuii'i' back- 
ground, which includes off-shell W*Z* 
and W*'y* production, as well as other 



1000 



Figure 17: The gluino pair production cross section, 
along the good-relic-abundance slice at tan /3 = 50 and 
mu > 0, compared to the total neutralino-chargino pro- 
duction cross section. On the upper x-axis we also in- 
dicate the reference gluino mass. 



diagrams |l68|. We have used Isajet 7.72 
to calculate the tt, ZZ and WZ with 
Z ^ TT backgrounds. For off-shell 
W*Z* and W*Y background calcula- 
tion, we have employed an exact tree 
level evaluation of the 2 — > 4 processes 
using Madgraph|^^ at the parton level. 
A similar calculation has been performed 
in Refs. pH] , where soft lepton px cuts 



were invoked [ 70 1 to try to maximize the 
signal when leptons originate from r 
decays. In our case, in the LHC envi- 
ronment, we will require three isolated 
leptons each with > 10 GeV and \ri\ < 2.5 throughout our analysis. In Ref. [^], also 

i. an invariant mass cut of m{l'i') < 81 GeV was invoked to reduce BG from the on-shell 
Z boson contribution, ii. m{i'i') > 20 GeV was used to reduce BG from the photon pole, 
and in. a transverse mass veto of 65 GeV < Mt{(.,i'i) < 80 GeV was used to reduce BG 
from on-shell W contributions. Finally, > 25 GeV was required. These cuts, dubbed 
SC2 by the Tevatron Run2 study group, will be invoked here, along with the somewhat 



stronger lepton px cuts. The BG rates after the SC2 cuts are listed in Table [4.1.2| . The 5a 
signal for cuts SC2 assuming 100 fb~^ of integrated luminosity is thus 1.38 fb. In Table 
4.1.2| , we also list a signal point with Mi = 110 GeV, where gluino pair production is still 
large and the signal is very robust. 

In Fig. 18, upper frame, we show the total pp — > W1Z2 cross section, and in the 
lower frame, the clean trilepton cross section after cuts SC2 along the line of ^^^^^ = 0.11 
as a function of Mi. Also shown by the horizontal mark is the 5a limit for 100 fb~^ of 
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Figure 18: (Upper) Total cross section for pp SUSY particles at the CERN LHC along the 
n^^/i^ = 0.11 line, as a function of Mi. (Lower): LHC clean trilepton cross section after cuts SC2 
and corresponding 5a discovery limits for 100 fb~^ integrated luminosity along the ^^_^h? = 0.11 
line as a function of Af i . 

integrated luminosity. We see that the CERN LHC reach for clean trileptons is possible 
out to Ml ~ 220 GeV, which is, in fact, comparable to the reach via conventional cascade 
decay signatures shown in Fig. left. 

While the results of Fig. |l8|, left, are valid along the line ^^^h^ = 0.11, it is also 
possible that ^^^h?' < 0.11 in scenarios with mixed dark matter. In this case, values of ^ 
smaller than those used in Fig. left, are possible. As becomes smaller, then the Zi 
becomes even more higgsino-like, and the relic density drops. The masses m^^ ^ and rriy^^ 
drop as well, as does the mass gap m^^ —m^^. The situation is illustrated in Fig. |l8|, right, 
for the case of Mi = 220 GeV, where various -ino masses are plotted versus a variable /x 
parameter. 
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Figure 19: (Left): Masses of lighter gauginos for a fixed Mi = 220 GeV along the region with 
^Zi^^ — O-ll- (Right): LHC total SUSY particle cross section (upper) and clean trilepton cross 
section after cuts SC2 (lower) along with corresponding 5cr discovery limits for 100 fb^^ integrated 
luminosity for a fixed Mi = 220 GeV along the line of varying fi. 



Thus, as n drops to smaller values, many of the -ino production cross sections rise. 
However, the trilepton energy and momentum distributions will diminish, in part due to 
the reduced parent particle masses, and in part due to the reduced sparticle decay mass 
gaps, which lead to smaller energy release in the chargino and neutralino decays. Thus, 
as the value of /i is reduced, production cross sections increase, while detection efficiency 
decreases. The trilepton cross section after cuts SC2 is shown versus /i for fixed Mi = 220 
GeV in Fig. |l^, right. Here we see that at large ^ values, the trilepton cross section 
after cuts is at the edge of observability. However, as /i decreases, the reduced detection 
efficiency wins out over the increasing production cross section (shown in the lower frame of 
Fig. ^), left, resulting in an overall diminished trilepton cross section. Thus, the trilepton 
cross section is actually maximal along the line ^^_^h'^ = 0.11, and diminishes in regions 
where O.^ < 0.11. 

In Fig. |2^, we show the opposite-sign/same fiavor dilepton invariant mass distribution 
from clean trilepton events at the CERN LHC for the case where Mi = 110 GeV. It is 
important to note two distinct mass edges in the plot. The first occurs from the kinematical 
edge from Z2 Ziff decay, and occurs at m^^ — m^^ = 56.4 GeV. The second comes 
from Z3 Ziff decay, and occurs at m^^ ~''^Zi ~ ^'^'^ GeV. The latter mass gap is close 
enough to the Z-pole that decay matrix element effect skews the invariant mass towards 
the high end of the range. 
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4.2 Prospects for sparticle detec- 
tion at the ILC 

The proposed International Linear Col- 
hder is projected to operate initially at 
y/s ~ 0.5 TeV with an integrated lumi- 
nosity of ~ 100 fb^^. In its later stages, 
the CM energy should be upgraded to 
~ 1 TeV. In Ref it has been 

shown that the reach of a linear collider 
can exceed that of the CERN LHC in 
the HB/FP region. This is because 
is always small in the HB/FP region, 
which forces detectable charginos and 
neutralinos- which should be readily ac- 
cessible to e^e~ colliders if the beam en- 
ergy is sufficiently high- to be relatively 
light, even if the scalars and the gluino 
are relatively heavy. 



Figure 20: Distribution in opposite sign/same flavor 
dilepton invariant mass from clean trilepton events at 
the CERN LHC for Mi = 110 GeV along the line of 
0.11. 



^1 



In Fig. |21|, we show total produc- 
tion cross sections for various -ino pair 
production reactions versus Mi along 

the line of constant O;; /i^ = 0.11. The 

^1 

left-hand frame is for a ^/s = 0.5 TeV 
collider, while the right-hand frame is for a 1 TeV collider. By examining frame a), we see 
that for Ml ~ 100 — 200 GeV, the -ino production reactions are dominated by chargino 
pair production, although a variety of other reactions such as Z1Z3 and W1W2 may 

also be accessible. The mass spectrum shown earlier in Fig. ^ shows that rn;^^ and 
are of order M2, and so should be heavy enough that two-body decays are accessible. How- 
ever, Wi, Z2 and Z3 have large higgsino components and are relatively light; they should 
typically decay via three-body modes, where the branching fractions are dominated by the 
W or Z boson propagators (since scalars are assumed quite heavy). Thus, the decays of 
Z2 and Z^ should be very similar in that Z2 — > Ziff and Z3 Ziff, aside from the size 
of the Z^ — Zi vs. the Z2 — Zi mass gaps. We see that the ultimate reach of the ^/s = 0.5 
TeV machine is determined by the ZiZ^ cross section. Z1Z2 production is more favorable 
kinematically, but has a lower total cross section due to a suppressed ZZ1Z2 coupling. The 
reach of a 0.5 TeV ILC along the ^^_^h'^ = 0.11 line is out to Mi ~ 220 GeV (corresponding 
to rrig ~ 1650 GeV), and is slightly below the reach of the CERN LHC. 

In frame b), for a ^/s = 1 TeV machine, we again see that W^Wi production is 
dominant over most of the range of Mi , while again a variety of other -ino pair production 
reactions should in general be present. In this case, the ultimate reach is determined by the 
Z1Z2 production reaction, and extends out to Mi ~ 500 GeV (corresponding to rrig ~ 3.5 
TeV), far beyond the reach of the CERN LHC. 
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4.3 Gluino lifetime in the HB/FP region 

The gluino decay width for g qqZi can be expressed as 

~ c — 1, 4.6) 

487r m~ 

where c represents a suitable combination of the neutralino-squark-quark couplings. The 
gluino lifetime, taking into account that the factor ~ 6.6 x 10~^, evaluating and a 
at Mz, can thus be cast as 



lQ-19 m| 

Conservatively assuming that c < 0.1, we obtain 



r ~ 1 sec (4.7) 



4 

r < 10"^^— I sec (4.8) 

In the case of minimal supergravity, we can draw the following general upper limit on the 
gluino lifetime. In mSUGRA, we always have rrig > 2.5 x mi/2, and, in the focus point 
region, where the gluino lifetime is maximal, we can safely take rriq ~ mg < 100 x 
We therefore have 

r < 10-" ( ^) sec. (4.9) 

Since in the focus point region mi/2 ^ 100 GeV, we find that, in minimal supergravity, 
r > 10~^^ sec. In order to detect a displaced vertex, the lifetime of a quasi-stable particle 
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Figure 22: The gluino lifetime as a function of the GUT-scale universal scalar mass mp, along 
mSUGRA slices at tan/? = 50, n > 0, Aq — for two fixed values of mi/2 = 500 and 1000 GeV. 



should be at least larger than rjisp ~ 10^ sec. This therefore entails that in mSUGRA 
gluinos are never "stable" inside a detector, which means that if a "stable" gluino , or a 
displaced gluino vertex is detected, the underlying SUSY theory cannot be mSUGRA. 

The largest gluino lifetimes are obtained in the focus point region at low values of 
mii2 and of tan/?, using a large top mass input. The spread in the gluino lifetime within 
mSUGRA can be significant, ranging from a minimum in the stau coannihilation region, 
where the squarks are the lightest possible, to a maximum in the HB/FP region. 

In the stau coannihilation region, we numerically find that (? ~ 0.5, which gives an 
estimated gluino lifetime between 10~^^ and 10^^^ sec, depending on the value of tan/3. 
The same range of gluino lifetimes is expected in the bulk and funnel regions. In the 
HB/FP region, one instead always obtains gluino lifetimes larger than 10~^^ sec. Since 
the gluino hadronizes if r > Thad ~ 6.6 x 10~^^ sec/(m7r/GeV) ~ 10~^^ sec, this means 
that the cosmologically viable parameter space of mSUGRA is split, from the point of view 
of gluino lifetimes, in a non-hadronizing gluino branch (coannihilation, funnel and bulk 
regions) and in a hadronizing gluino branch (HB/FP region). A hadronizing gluino gives 
rise to i). additional hadrons through fragmentation, ii). smeared out spin correlations 
and Hi), reduced energy in daughter particle energy distributions. All of these effects 
will, however, be very difficult to determine experimentally in the environment of an LHC 
detector. 
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5. Conclusions 



The HB/FP region of SUGRA models is very compelling from both a theoretical and 
phenomenological viewpoint. The multi-TeV scalars act to suppress unwanted FCNCs, 
CP-violating dipole moments and proton decay, while maintaining low fine-tuning. In ad- 
dition, the HB/FP region is entirely consistent with the WMAP CDM constraint, due to 
the presence of a mixed bino-higgsino LSP. Previously, many investigations of phenomenol- 
ogy in the HB/FP region took place in the mSUGRA (ttt-q, 'fn^i'i) plane. Calculations of 
HB /FP phenomenology in terms of these GUT scale parameters are plagued by numerical 
and parametric uncertainties. We propose here, instead, a presentation of HB/FP phe- 
nomenology in terms of weak scale (Mi, [i] parameters, wherein i). the entire HB/FP 
region can be displayed, including the — > region, and O). the parameter plane is much 
more stable against variations in other parameters such as rrit and tan /?. 

We present the HB/FP in the (Mi, [i] plane, and show regions that are allowed and 
disallowed by theory, collider and astrophysical constraints. We also show prospects for 
future DM searches in the (Mi, n) plane. Prospects for exploration of the WMAP allowed 
portion of the HB/FP plane are excellent for Stage 3 direct dark matter detectors, for 
indirect detection via neutrino telescopes, and for antideuteron searches via the GAPS 
experiment. Prospects may also be good for 7 ray searches from neutralino annihilation in 
the galactic core, but these estimates depend sensitively on the assumed dark matter halo 
distribution near the galactic center. 

Regarding collider searches in the WMAP allowed portion of the HB/FP region, we 
find that the LHC, with 100 fb~^ of data, should be able to probe to ~ 1.8 TeV 
via conventional gluino cascade decay signatures. We also studied the reach of the LHC 
for SUSY in the clean trilepton channel, including backgrounds from W*Z* and W*j* 
production. In this case, the LHC reach for 100 fb^^ is slightly smaller, out to rrig ~ 1.65 
TeV. A linear e^e~ collider will have a reach in the HB/FP region determined by Z1Z2 
and production. A y/s = 0.5 TeV LC will have a reach to rUg ~ 1.65 TeV, while a 

-y/s = 1 TeV LC will have a reach beyond that of the CERN LHC, out to rrig ~ 3.5 TeV. 
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